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Quantum Mechanics as Asymptotics of 
Solutions of Generalized Kramers Equation 

q ■ E. M. Beniaminov 

o 

oo 

r*j ' We consider the process of diffusion scattering of a wave function 

Oh given on the phase space. In this process the heat diffusion is con- 

^ (— | sidered only along momenta. We write down the modified Kramers 

equation describing this situation. In this model, the usual quantum 
description arises as asymptotics of this process for large values of re- 
sistance of the medium per unit of mass of particle. It is shown that in 
this case the process passes several stages. During the first short stage, 

■ the wave function goes to one of "stationary" values. At the second 
long stage, the wave function varies in the subspace of "stationary" 
states according to the Schrodinger equation. Further, dissipation of 

I/"") ■ the process leads to decoherence, and any superposition of states goes 

■ to one of eigenstates of the Hamilton operator. At the last stage, the 
mixed state of heat equilibrium (the Gibbs state) arises due to the 
heat influence of the medium and the random transitions among the 
eigenstates of the Hamilton operator. 

Besides that, it is shown that, on the contrary, if the resistance of 
^ ' the medium per unit of mass of particle is small, then in the considered 

■ model, the density of distribution of probability p = \tp\ 2 satisfies the 
standard Liouville equation, as in classical statistical mechanics. 



1 Description and some properties of the 
model 

As in [2], we consider certain mathematical model of a process whose state 
at each moment of time is given by a wave function, which is a complex 
valued function <p(x,p), where (x,p) G R 2n , on the phase space, and n is 
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the dimension of the configuration space. In contrast to quantum mechanics, 
where the wave function depends only on coordinates or only on momenta, in 
our case the wave function depends both on coordinates and on momenta. By 
analogy with quantum mechanics, it is assumed that wave functions obey the 
superposition principle, and the density of probability pu(x,p) in a bounded 
domain of the phase space (x,p) G D C R 2n , corresponding to the wave 
function tp(x,p), is given by the standard formula 

Pd{x,p) = \(p(x,p)\ 2 / / \Lp(x,p)\ 2 dxdp. (1) 

In quantum mechanics, the time evolution of the wave function can be 
defined by the Feynman path integral [3j. The Feynman principle assumes 
that if at the initial moment t = to a wave function ip(xo,Po, to) is given, then 
the value of the wave function at the point (x,p) at the moment t = t\ is 
defined by the integral over all paths {x(t),p(t)} joining the points (x ,Po, to) 
and (x,p,ti), of the quantity exp (— | [V(x(t)) — p 2 (t) / (2m)]dt^j , where % 
is the Planck constant, with respect to certain "measure" on paths defined 
by Feynman. 

In contrast to Feynman's assumption, in the present paper we study the 
model in which the Feynman measure on paths is replaced by the probability 
measure of the diffusion process (the heat Brownian motion) given by the 
Kramers equation [1], [5]: 

dl^y(9Vdl_P 1 dl) +1 y±( p . f + kTm K) , (2) 

dt ~[\dxj dpj mdxjj J^idpjK 3 dpj)' 

where f(x,p,t) is the density of probability distribution of the particle in 
the phase space at the moment of time t; m is the mass of the particle; 
V(x) is the potential function of the external forces acting on the particle; 
7 = (3/m is the resistance coefficient of the medium per unit of mass; k is 
the Boltzmann constant; T is the temperature of the medium. 

This is the classical Kramers equation describing the diffusion motion 
of a particle in the phase space under action of external forces defined by 
the potential function V(x), the heat medium with temperature T, and the 
medium resistance per unit of mass 7. 

Consider the following modified Kramers equation for the wave function 
cp(x,p,t): 

^ = A V + 7 I^, (3) 
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" (dV dip p,dcp\ if 2 „ Api\ ... 
where Ap = " 5 (™? + V - g ^ (4) 

( 5) 

Equation ([3]) is obtained from the Kramers equation (J2]) by adding to the 
right hand side the summand of the form —i/h(mc 2 + V — p 2 /(2m))(p and 
by replacing multiplication of the function ip by pj in the diffusion operator 
by the action of the operator (pj + ihd/dxj) on the function ip. 

Adding the summand —i/h(mc 2 + V — p 2 /(2m))p is related with the 
additional physical requirement that the wave function at the point (x,p) 
oscillates harmonically with the frequency l/h(mc 2 + V — p 2 /(2m)) in time. 

The requirement of harmonic oscillation of the wave function <p at the 
point (x,p) with large frequency l/h(mc 2 + V — p 2 /(2m)), in the case when 
mc 2 is much greater than V, gives that the shift of the wave function with 
respect to coordinate Xj with conservation of the proper time at the point 
(x,p) yields the phase shift in the oscillation of the function if. This also 
gives that the operator of infinitesimal shift d/dxj is replaced by the operator 
d/dxj —ipj/h. (For a more detailed explanation, see 0.) Respectively, if one 
multiplies this operator by ih, then one obtains the operator pj + ihd/dxj 
used in the modified diffusion operator B. 

Let us proceed to the study of equation ()3]). 

In order to separate mathematics from physics in this equation, let us 
make a change of variables and let us pass to dimensionless quantities: 



/ / V i VkTm Tr ,. , V(x) 

In the new variables equation ([3]) takes the following form: 

dip kT . . 

^ = ^V + ^V, (7) 
where A ip = ^( t^-ttt ~p[ 



■J x V dx'j dp'j 3 dx'j 
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The parameter of the model described by equation (j7j), is the dimension- 
less quantity kT/(jh), which we denote by s. 

Let us assume that e = kT/(jh) is a small quantity which is the small 
perturbation parameter in equation © with the non-perturbed equation 
dip/dt' = B'ip, i. e., the equation 




Note that the smallness of the quantity e = kT / ( r yh) requires that the friction 
coefficient of the medium per unit of mass 7 = (3/m = (k/h)T/e = 1.3 • 
10 u T/e be greater than 1.3 • 10 n T, since k/% = 1.3 • 10 11 . 

Let us substitute into equation (|TU|) the Fourier integral presentation of 
<p(x',p',t') with respect to x': 

<P&>M = J R J(s',p',t')e^'ds', (11) 

where (p{s' ,p' ,t') = — ^ <ptf tI /, t')e- is ' x ' dx' . (12) 
We obtain that (p(s',p',f) satisfies the equation 




The operator of the right hand side of this equation is well known (see, 
for example, [6]). This operator has a full set of eigenfunctions in the space 
of functions tending to zero as \p'\ tends to infinity. The eigenvalues of this 
operator are nonpositive integers. The eigenfunctions corresponding to the 
eigenvalue have the form 

~ / / a 1 7/ a - (p'-"') 2 

^(s,p) = j^y^Hs)e * , 

where ip(s') is an arbitrary complex valued function of s' G R n . 

The rest of eigenfunctions are obtained as derivatives of the functions 
L Po( s ',p') with respect to p', and have eigenvalues —1,-2,..., respectively, 
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depending on the degree of derivative, and the projector P to the subspace 
of eigenfunctions with eigenvalue has the form 

1 ~ <p'-s') 2 ~ r 

<fo(s',p') = P <p = p3^ ^(s')e ~ , where i[)(s') = <f(s',p')dp'. 

(14) 

Hence, considering equation (fl3l) in the basis of these eigenfunctions, we ob- 
tain that each solution <p(s', p', t') of this equation tends exponentially in time 
with exponent —1 to a stationary solution of the form (p . Therefore, taking 
into account the presentation ffTTj) of the function (p(x',p',t) via (p{s' ,p' ,t'), 
we obtain that "stationary" solutions <f (x',p') of equation ([7]) look as follows: 

Let us present the function ip(s'), in its turn, as the Fourier integral: 



Substituting this presentation into the preceding expression and integrating 
over s', we obtain 

Mx',p) = j^j, J R2 J(y')e-^e^'^ds' dy> 
1 ' ^(y')e-^^'-yyy' 



(2ir) n Jr" 
or, taking into account f|T4|) . 



= Po<P = ^-^{y')e- {d ^^^-^dy\ where ^(y') = J y{y',p')dp'- 

R n R n 

(15) 

Thus, if e is small, then at the time i! of order 1, a solution of equation ([7]), 
starting from an arbitrary function (p, will become close to a function of the 
form ifQ which, in the initial coordinates ([6]), reads: 



, , kTm(x-yf\ ( ip(x - y)\ 

Mx ' p) = (2^f L m exp ( — w — J exp {—fr- 1 ^ 
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Later the solution of equation Q, evolving in time, will be close to the 
subspace of "stationary" functions of the form ip . 

Note that the "stationary" functions obtained here coincide, up to a nor- 
malizing factor, with the "stationary" functions of Theorem 1 obtained in 
the work |2], if one puts kTm/Ti equal to b/a, where a 2 and b 2 are the dif- 
fusion coefficients with respect to coordinates and momenta in the model 
considered in [2J. Hence some results obtained in [TJ, [2J, hold true also in 
our case. We shall cite them without proof. 

The results are obtained by the perturbation theory method up to second 
order, of equation (JTj) dtp/dt = sA'ip + B'lp for e < 1, where A' is a skew 
Hermitian operator and B' is an operator with nonpositive discrete spectrum. 

Note that this equation does not preserve the norm of the function ip, 
which we normalize to obtain the distribution p(x,p). The corresponding 
equation preserving the norm, which we actually study, reads dip/dt = eA'tp+ 
B'ip + kip, where k = — ((B'p, ip) + (ip, B'<p))/(2(<p, ip)), but it is nonlinear. 

2 The main results 

Denote by H(x,p) = p 2 /(2m) + V(x) + mc 2 the Hamilton function of the 



Theorem 1. The motion described by equation (T3)] ; for small e = 
kT / (jh) asymptotically splits into rapid and slow motion. 

1) After rapid motion the arbitrary wave function ip(x,p,0) goes, in time 
of order 1/7, to a function which after normalization has the form 



where M\ = 1 and X (x - y) = e -kT m { x - y ?/{2tf) ^ (17) 



The wave functions of the form j[Tb)) form a linear space. The elements of 
this subspace are parameterized by the wave functions if>(y) depending only 
on coordinates y e R n . 

2) The slow motion starting from the wave function (po(x,p) of the form (To)) 
with nonzero function ip{y), goes inside the subspace and is parameterized by 
the wave function ift(y,t) depending on time. The function ip(y,t) satisfies 



system. 




(16) 
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the Schrodinger equation of the form ihdip/dt = Hip, where 



h 2 / n f) 2 ih\ kT 
H * = ~ (E ^r) + V(y)il> + \ni> + mcty (18) 

Proof of the first part of Theorem 1 is given in [2] . Proof of the second 
part of the Theorem is given in Appendix 1 to the present paper. 

Theorem 2. The projection operator P transforming an arbitrary in- 
tegrable function ip(x,p) on the phase space into the function of the form 
(To)) (but without normalization), obtained after rapid motion described in 
Theorem 1, reads as follows: 

1 r _kT(x-x') 2 ip(x-x') r 

PoV = (2 h)n J if>{y)e ^ e » dy, where ip(y) = J <p(y,p)dp. (19) 

^ ' R n R 

This Theorem follows from the definition of the operator Pq and from 
formula (115j) expressed in the initial coordinates fl6]). 

Theorem 3. If ip(x) is a wave function on the configuration space and 
<Po(x,p) is the wave function on the phase space corresponding to it by for- 
mula l[Tb)) . then the density of probability p(x,p) = \<po(x,p)\ 2 in the phase 
space is given by the following formula: 



^• p) =(2^u^j r\ x+ — r\ x+ — ) (20) 

/ kTm(x") 2 \ ( kTm(x') 2 \ (ix'p\ , .. , , 
exp ^— exp ^j- exp — \dxdx. 



In contrast to quasidistributions 

W(x, p) = j^j* (* ~ ~) r (x + | J exp j dx' 

defined by Wigner [7], the density p(x,p) in the phase space, given by the 
expression (1201) . is always nonnegative. Its expression differs from the ex- 
pression of the Wigner function by exponents under the integral, which yield 
smoothing with respect to distribution densities close to the delta-functions. 
Proof of Theorem 3 is given in [2] . 
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The algebra of observables given by real functions on the phase space, 
averaged by densities of probability distributions of the form (I2"0"j) . has been 
studied in [8]. 

Since the probability distribution p(x) in the configuration space is ex- 
pressed by the formula p(x) = J R3 p(x, p)dp, by integrating the expression (120]) 
over p we obtain the following statement. 

Corollary 1. If ip(x) is the wave function on the configuration space, 
then the corresponding density p(x) in the configuration space is given by the 
formula 

p(x) = f My)\ 2 x 2 (x-y)dy, (21) 

where x( x ,y) expressed by relation fT7| ]. That is, p(x) is obtained from 
\ip(x)\ 2 by smoothing (convolution) with respect to the density of the normal 
distribution with dispersion h 2 / (2kTm) , and the exactness of defining coor- 
dinate is bounded by the quantity ~ Ti/ \/2kTm called the de Broglie length of 
the heat wave. 

Theorem 4. // the number e = kT/^h) <C 1, then the first order per- 
turbations of the zero eigenvalue of the operator B' in equation equal to 
the eigenvalues of the operator —i/{^fh)H, where H is the Hamilton oper- 
ator given by formula If the real parts of the second order perturba- 
tions corresponding to these first order perturbations are different from one 
another, then any solution of equation (Gj) will go in time proportional to 
l/(7£ 2 ) = 7^ 2 /(fcT) 2 ; to one of the eigenstates of the Hamilton operator. 

Proof of the Theorem is given in Appendix 2. 

This theorem describes the process called in the literature by decoherence 
of quantum states [9], [10], [TTJ. The form of the estimate of the time of deco- 
herence given in Theorem 4, somewhat differs from the form of the estimate 
given in the literature, for instance in [10J. The study of the correspondence 
of these estimates is the subject of another future work. 

According to Theorem 4 the process described by equation ([3]), in time 
proportional to , yh 2 /(kT) 2 , goes to one of eigenstates of the Hamilton (en- 
ergy) operator. Further, on large scales of time the system, under the action 
of the heat medium, will jump from one eigenstates to others due to large 
deviations of the random process. In the limit as t — > oo the system will go 
to the mixed state corresponding to the heat equilibrium Gibbs state. 

We have studied equation ([3]) for small value of e = kT / (•yh) and obtained 
that the process described by this equation is asymptotically close to the 
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process having the standard quantum description. Let us now consider the 
same equation in the case when the quantity e is large. 

Theorem 5. If e = kT/(jh) > 1, i. e. jh/(kT) < 1, where 7 = /3/m 
is the resistance of medium per unit of mass, then the operator B in equa- 
tion (T3|) can be neglected, and in this case the density of probability distribu- 
tion p(x,p,t) = <f(x,p,t)(p*(x,p,t) satisfies the following classical Liouville 
equation: 



Proof of the Theorem follows from the fact that equation ([3]) without the 
operator B is a partial differential equation of first order, consisting of sum 
of the Liouville operator and the operator of multiplication by a function. 
The solution (p(x,p,t) of such equation is obtained from the initial state 
<p(x,p) = <p(x,p,0) and the characteristics x(t),p(t) of the equation with the 
following initial conditions: x(0) = x, p(0) = p, in the following form: 



p(x,p,t) = <p(x,p,t)<p*(x,p,t) = <p(x(t),p(t))<p*(x(t),p(t)) = p(x(t),p(t)). 



Therefore, the phase of the wave is inessential, and the density p(x,p,t) 
satisfies the Liouville equation. 

Thus, Theorem 5 states that for small value of (3h/(kTm) the process 
described by equation (j3J) is asymptotically close to the process with the 
classical (non-quantum) description of the motion of the particle. This case 
arises when the mass m of the particle is large relative to the medium resis- 
tance (3. 

3 Conclusion 

In quantum optics, one now widely uses methods of study of dynamics of 
quantum processes in the phase space, see, for example, [12J, where, in par- 
ticular, one considers dynamics of wave packets and interference in the phase 
space. On this way, one manages to model dynamics of ions in traps, optics 
of atoms in quantum light fields, etc. "This approach which stresses the 



dp 
di 




(22) 



<p(x,p,t) = tp{x{t),p{t))exp (~i/hj (H(x(t),p(t))-p 2 (t)/(2m))dt 



Respectively, we have: 
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fundamental role of phase variables allows one to expose and interprete very 
clearly various branches of quantum optics..." (from the abstract of the book 

H21). 

These achievements lead to an assumption that it would be useful to 
consider not only the behavior of distributions in the phase space, as it is 
done, for example, in [T2], [13], [33], but also to introduce the wave function 
itself on the phase space. 

In the present work, following this direction, we construct a diffusion 
equation for the wave function with large frequency of oscillations in the 
phase space, describing the process of heat scattering of a wave in the phase 
space. In the presented model one meets both classical and quantum me- 
chanics behavior of the particle. If the quantity e = kTm/(/3h) is small, then 
in this model the behavior of the particle, after short transition stage (of 
order m/(3 < Ti/kT = 0.77 • 10~ 12 /T sec.) amounts to the standard picture of 
quantum mechanics with the Heisenberg indeterminacy principle and with 
the Schrodinger equation describing the dynamics. And if the quantity e 
is large, then the particle behaves according to classical mechanics, and the 
density of probability distribution of the particle in the phase space satisfies 
the classical Liouville equation. 

In the general case, the behavior of the particle described by this model 
is of mixed nature. It would be interesting to analyze the model in this case 
and compare it with results of experiments for particles with intermediate 
values of e. 

One should also compare with experimental data the estimate of time of 
decoherence given in Theorem 4. Besides that, one should find a theoretical 
estimate of the transition time of the quantum system to the mixed state of 
heat equilibrium in this model. 

One should also acknowledge numerous works of predecessors, due to 
which the subject of this paper arose, the problem has been stated and the 
methods of its solution have been found. This is a separate large work. 
In science, posing right questions gives no less than the results obtained. 
Bright examples of it are the questions of A. Einstein, due to which quantum 
mechanics was founded and develops up to now. 

In the works of V. P. Maslov [151 EE] one has already studied certain 
problem of description of motion of a distribution of charges in the phase 
space under action of a random field. In [15] under certain assumptions it is 
proven that if the initial distribution of charges, depending on coordinates 
and momenta, belongs to certain subspace, parameterized by complex valued 
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functions depending on coordinates, then the distribution does not leave this 
subspace, and the dynamics of such system is described by the corresponding 
Schrodinger equation. This result has certainly influenced the author while 
posing the problem of the present work. 

Acknowledgements: I am deeply grateful in my heart to Professor 
G. L. Litvinov for the many year attention to my work, for understanding, 
and for help in the right organization and exposition of the material. I am 
grateful to professor A. V. Stoyanovsky, who translated this paper to English. 
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Appendices 



Appendix 1. Proof of Part 2 of Theorem 1 

Consider equation ([3]) in the dimensionless system of variables (Q. In these 
variables, the equation takes the form ([7]), and "stationary" solutions, to 
which arbitrary solutions of equation ([7]) tend at time t' of order 1, have the 
form ( 1T5|) . 

Let (p (x',p', t') be a function of the form ([15]) corresponding to the func- 
tion i])(y\t'). Let us substitute this expression into equation ([7]), and let 
us take projection of both parts of this equation to the space of functions 
ip(y' ',t') by formula (|T5|) . We have: 



or, taking into account that B'tp = 0, after substitution of expression 
<po{%',p', f) in the form ( fT5l) . we obtain: 



(2ttW dt 



ydp 

kT 1 



7^ (2tt) 



A'^{y', t')e—^e iv ' {x '- vl) dy'd P '. 



Let us integrate the right hand side of this equality over p' and over y'. 
Noting that in the right hand side of the equality we have the delta function, 
we obtain: 

\ = — - — — A'ib(y',t')e 2 e w {x y) dy'dp'. 

V ^ R2n 

Taking into account expression (jBJ) for operator A', we deduce from the latter 
equality that 

foj) = KT 1 rf»/dV'd ,d\ .fmc 2 " jp'j) 2 

m'~ ^h{2n) n J 2 \p l \dx' j dp' j Pj dx' j ) \kT + p x 2 



kT 



x^(y',t')e-'-^e ip '^~y^dy'dp' = —{h + J 2 + J 3 + J 4 ), (23) 
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where 



(2ii) n J £rl Bx'a dp 

1 r " , d ( . , , ,, s _ ± 



nj (^(y'.Oe-^c*'^^) W; (25) 

+ V\x')^{y\ty- {J ^e ip ' {x '- y ' ) dy'dp'- 1 (26) 



1 f ( mc 



(27r) n J V kT 



i 1 



,,,)■>/ E^^d/'^Oe-^e-'^'-^W. (27) 

Consider the integral I\ given by expression (|24p . Let us exchange sum- 
mation and integration, let us put behind the sign of integral expressions 
not depending on integration variables, let us compute the derivatives with 
respect to p'j, and let us integrate the remaining integrals over p' and y'. We 
obtain: 



^TtV^^T 1 / W,ty-^(x' 3 -y'^'-y'Uy'dp' = (28) 
l Z7r J j= i i ^ 2n 



Consider the integral I2 given by expression (125 p . Let us exchange sum- 
mation and integration, let us transfer the derivatives with respect to x'j 
behind the sign of integral, let us replace the expressions p'j exp(ip'(x' — y')) 
by equal expressions id exp(ip'(x' — y')) / (dy'j) , and let us integrate the ob- 
tained integrals by parts. We have: 

1_ f-JL [ i>(y>,t>)e-^iJL e i P '(*'-v') dy > dp > 



E^7 / v(i/ / ,t / )(4-^) e_ ^ eV(a: '^ ) W 
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Consider the integral /3 given by expression ( I26p . Let us transfer behind 
the sign of integral the expressions not depending on the integration variables, 
and let us integrate the remaining integral over p' and y' . We obtain: 



—i 



-<(^ + VV))^(^0- (30) 



Consider the integral I4 given by expression (1271) . Let us exchange sum- 
mation and integration, let us transfer 1/2 behind the sign of integral, let 
us replace the expression (p'j) 2 exp(ip'(x' — y')) by the second derivative of 
the function — exp(ip'(x' — y')) with respect to y'p and let us integrate the 
obtained integrals by parts. We have: 



n r ,_ , 2 pp. 

j^E / WJ)e-^^e^'-y'Wdp> 

,E/(^-2^-^) + ^-^ + ; 



xe-^e^'-y'^dy'dp' 

Let us substitute the obtained expressions for integrals J 4 into 

equality (12*3"]). let us sum up similar terms, and let us transfer —i behind the 
brackets. We obtain, 

dtb i kT ( 1 d 2 Ir . mc 2 n\ , . , .. 

If in the obtained equality we pass to the initial coordinate system ([6]), 
then we obtain the equality f fl8|) required in Part 2 of Theorem 1. 
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Appendix 2. Proof of Theorem 4 



Consider equation (J7J) of the form dip/dt' = eA'ip + B'ip. Let P Q be the 
projector onto the subspace of eigenfunctions of operator B' with eigenvalue 
0. Then, by definition of Pq, we have the equalities: 

P P = P ; PoB' = B'P = 0. 

We shall denote by <po a function belonging to the image of projector P , i. e. 
V? = Po<Po- 

Note that by construction of operator H from Theorem 1, given in Ap- 
pendix 1, the operator ePqA' on the subspace of functions (fo, in the presenta- 
tion by the functions ip{y) has the form —ij ( r yh)H . Hence the eigenvalues of 
operators ePqA' and —i/ (jf^H coincide. Since the operator H is self-adjoint, 
it has a complete system of eigenfunctions, and therefore the operator ePqA' 
on the subspace of values of the projector P also has a complete system of 
eigenfunctions. 

Let us state the result thus obtained Lemma. 

Lemma. The eigenvalues of operators ePqA' and —i/( , ~fh)H coincide. 
The operator ePqA' has a complete system of eigenfunctions on the subspace 
of values of the projector P . 

Consider the eigenvalue problem for equation (j7j) of the form \ £ {p £ = 
(eA' + B')<p £ . 

According to perturbation theory, let us look for solutions as series: 

X £ = A + eAi + e 2 \ 2 + ... 

tp £ = ip + ELfi + E 2 Lf 2 + ••• 

Let us substitute these series into the eigenvalue equation and compare co- 
efficients before equal powers of e. We obtain: 

(e ) -Wo = B'tpo; 

(ei) Ai^o + Ao<^i = A'(p + B'cp 1 ; 

We are interested in perturbations of eigenvalue Ao = 0. In this case 
equation (eg) means that (fo is an eigenfunction of operator B' with eigen- 
value 0. 

Let us apply the operator ePq to both parts of equality (ei). Taking into 
account the equalities PoV^o — V^O; A = 0, and PqB' = 0, we obtain: 

eX^q = eP A'(p . 
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Hence the quantity eAi is an eigenvalue of the operator ePqA' on the 
subspace of values of the projector P , and by the Lemma the quantity e\i 
is an eigenvalue of the operator —i/(j?i)H. On the other hand, the quantity 
e\\ is by definition a first order correction to the eigenvalue of the operator 
eA' + B'. 

Thus, the first part of Theorem 4 is proved. 

To prove the second part, let us present an arbitrary function <po(t') from 
the subspace of values of the projector P as a sum (or integral for continuous 
spectrum): 

oo 
k=l 

where <^o are the eigenfunctions of the operator PqA' with eigenvalues A^. 
This is possible, since by the Lemma such functions form a complete system. 
Besides that, the same Lemma implies that the numbers A^ are pure imag- 
inary, since they are eigenvalues of the operator —i/(e , yh)H, where H is a 
self adjoint Hamilton operator with real spectrum. 

Let Ag « e\\ + e 2 \2, and (p* « (p* + e<p\ be the approximations 
of eigenvalues and eigenfunctions corresponding to eigenfunctions (p^. By 
assumptions of Theorem 4 the real parts Re\\ are different for various k. 

Consider equation (j7|) restricted to the subspace spanned by the basis 
vectors of the form p)\ with the function of the form 

oo 
k=l 

In this basis the equation splits (up to terms of order e 3 ) into a system of 
equations numbered by k of the form: 



dt> 



( £ A^ + e 2 A^)c fc (0 



whose solution reads as follows: Cfc(t') = Cfe(O) exp((eAf + e 2 \2)t'). Hence, the 
solution of equation (JTJ) is approximately presented in the following form: 



^(O = E^(0)exp((eAt + e 2 A^K)^ 

oo 

= E c*(0) exp((eAt + e 2 \ h 2 )t'M + ep>\ 



k=i 
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£ Cfc (0) exp((eA* + £ 2 A 2 fc )t>o + E c *(°) exp(( £ A* + £ 2 A*)t')^ 



k=l 



k=l 



(32) 



where tp (t') = £ c fe (0) exp((eA^ + s 2 \ k 2 )t')^ (33) 
fc=i 

oo 

Mt') = E c *(°) exp((eA* + e 2 A§)0^- (34) 



fc=i 



Since e is small, equality ( 132]) implies that the solution ip(t') has a small 
difference with the function (fo(t') given by expression fl33l) . 

Since eA^ is imaginary, each summand in the sum giving the function 
ipo(t') in expression ( 1331) . decreases in time i' = 7t proportionally to 



where -Re^Jj) is the real part of the number A|. This implies that after the 
time t ~ l/(7£ 2 ) this sum will be determined by the summand with the 
maximal number i?e(A|) among the summands with the valuable coefficients 
c fe (0). 

Thus, Theorem 4 is proved. 
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